Review of [theoretical] models 


“State of the art of the models and simulations of 
galaxies in groups and clusters” 


Goal: 

= Overview of current models of different types 
Scope of physical ingredients in each 
(Limitations and inherent assumptions) 
Comparison of numerical approaches, regimes of 
applicability, classes of objects, etc. 


= Extensive review of the scientific results of the topics 
of this workshop as addressed by models 
= Focused on details of feedback (tomorrow) 


Review of [theoretical] models 


Bruno Henriques 

"= SAMs 

Stephanie Tonnesen & Elke Roediger 
= |dealized & wind-tunnels 
Annalisa Pillepich 

= Cosmological hydro sims 

Dylan Nelson 

= ..in the context of Illustris / TNG 
Yannick Bahé 

= ... in the context of EAGLE 
Collective discussion 


Review of [theoretical] models 


= Bruno Henriques 
= SAMs 


Challenges in galaxyformation 
theory 


Bruno Henriques (Zwicky Fellow, ETH Zurich) 


Simon White, Peter Thomas, Simon Lilly 
Raul Angulo, Scott Clay, Benoit Fournier, Qi Guo, 
Gerard Lemson, Volker Springel, Rob Yates 


Basic principle of SAMs: Gravity 


Content of the Universe known to a 
few percent!!! 


Dark Matter 


We have a working model of the 
universe in which 85% of matter Met 


only interacts through gravity 


Before Planck 


Dark Matter 


After Planck 


Evolution of structures is “Tully known” 


masses, sizes, temperatures, positions and 
velocities of dark matter haloes are “known” 


f NE 
d o rm J 
- ASS a 
DA te g X 
v & 7 
e Ti 
ý 5 O 
À 9 Q e Q « 3 
» L Ti 
x >. 
| | "4 
2 Q Q 
Wye: ? 
S, 2^ | N/ 
3 
AE 
på - 
T ‘a 
IY 
Å 
d 


Accretion 


baryonic mass is given by cosmology: 
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Cooling 


What happens to the baryons? hot atmosphere vs rapid cooling 
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Henri Thomas i I 
Henriques & Thomas (2010), Henriques etal. (2013); Henriques et al. (2015). 
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tidal and ram- 
æ pressure stripig 


= © Environment removes the fuel for star formation in satellite galaxies. It is the 
predominant quenching mechanism in satellite galaxies 


Lack of Primordial Infall - 


= * Environment removes the fuel for star formation in satellite galaxies. It is the 
predominant quenching mechanism in satellite galaxies 


Hot Gas Striping - 
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tidal and ram- 
pressure stripig 


= * Environment removes the fuel for star formation in satellite galaxies. It is the 
predominant quenching mechanism in satellite galaxies 


Tidal disruption of cold gas and stars 
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=" tidal disruption 


= * Environment removes the fuel for star formation in satellite galaxies. It is the 
predominant quenching mechanism in satellite galaxies 


Tidal disruption of cold gas and stars 


= * Environment removes the fuel for star formation in satellite galaxies. It is the 
predominant quenching mechanism in satellite galaxies 


Positions of orphans $ 


= © Environment removes the fuel for star formation in satellite galaxies. It is the 
predominant quenching mechanism in satellite galaxies 
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Environment removes the fuel for star formation in satellite galaxies. It is the 
predominant quenching mechanism in satellite galaxies 


AGN and Environment Quenching 


higher quenched fractions for higher higher 


passive fraction vs stellar mass 
stellar masses and denser environments 


data from Wetzel et al. 2012 - massive galaxies quenched due to AGN 


- most low mass galaxies are star-forming centrals 


- 40% of low mass galaxies are satellites of which 
~50% are quenched (20% quenched low mass) 


(a(Maow/Mo )=[14.5,15.0] passive fraction vs halo mass 
109(Mroio/ Mo )=[14.0,14.5] 
— 109(Mpgig/Mo )=[13.5,14.0] 
— 109(Mnoip/Mo )=[13.0,13.5] 
109(Mroio/ Mo )=[12.5,13.0] 
log(Mnaio/Mo )=[12.0,12.5] 
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Review of [theoretical] models 


= Stephanie Tonnesen & Elke Roediger 
= Idealized & wind-tunnels 


Gas-related environment effects on cluster 
galaxies = ram pressure stripping = gas stripping = ICM 
stripping 


Cluster General idea: ICM head wind 


strips galaxy' s gas. 
Does it work? Yes! 


Pram ~ 9 XY gas (Gunn&Gott 72) 


Eon n Picm = Pism 


Gas stripping — the devils in the details 


Most simple idea: Object in a wind, head wind 


Shape factor: sphere / disk / inclination angle 
Variation of head wind / ram pressure during cluster passage 


For disk galaxies: 
huge dynamic range in gas T, p, spatial and time scales 


Proper SF and stellar feedback 


Turbulent ISM, ICM, bulk motions in ICM 


Magnetic fields in ICM, in ISM 


"Gas model" for ICM? Hydro/MHD/Kinetic code or ???, two- 
temperature? Mixing? 


Proper radiative transfer to heat and cool the cool gas embedded in 
hot ICM 


Proper radiative transfer to make mock observations 


Idealized simulations 


e Single/few galaxies 
e Non-cosmological 


e Choose the physics 
implemented, well 
controlled! 


e Manual initial conditions 
(DM, stellar, gas, ICM 
wind, ...) 


e Scan parameter space 


proj. gas density m 


s- 7 e Can be high resolution 
Roediger+ 2014 


Consensus reached: 


e Gas stripping works, 
truncated galactic gas, 
head-tail structure 


S g | Model (edge-on | 

e Gunn&Gott criterion / 
pressure comparison is 
decent prediction of 
TRUNCATION of galactic 
gas, almost independent 
of inclination 


e Actual gas removal takes 
time — still-bound gas in 
tail region, fall-back 


Takeda+84, Stevenson+99, 
Abadi*99, Schulz+01, Vollmer+ ..., Innsbruck group, 
Tonnesen- ..., Roediger+..., and more! 
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»roj. gas density 


Harder questions: impact on SF 
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Less massive disk 
More gas-rich 

Different resolution 
Different SF recipe 
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ratio star formtaion to vacuum star formation 
Ja e 


surrounding gas density [g/cm "] 


Fig. 5. Ratio of the star formation of the different simulations to the star 
formation integrated over 500 Myr for the isolated evolving galaxy. 


Harder questions: impact on SF 


Tonnesen Kapferer Roediger 
No Enhancement factor 2-10, Strong enhancement 
Enhancement. A lot of SF in the tail of gas that is stripped 
immediately afterwards 
AMR SPH AMR 


cells in disk^"5x10" 
Mass refine: 4900 M,,, 
cooling: 300 K 

Toe: 1.1X104 K 

Psp: 3.85 x 102 g cm? 


particles in disk: 2x10? cells in disk~5x10° 
particle Mass: 3.4x10* M.n Fully refine to 30 kpc 
cooling: 10^ K cooling: 10^ K 

Toe: 106 K Top: 1.5x104 K 

Psr: “7 x 1028 g cm? (?) psg: 3. x 107^ g cm? 


consider resolution and star an 
formation recipes 


Tonnesen Kapferer Roediger 
AMR SPH AMR 
cells in disk"5x10" particles in disk: 2x10? cells in disk~5x10° 
Mass refine: 4900 M,,, particle Mass: 3.4x10^ M,,, Fully refine to 30 kpc 
cooling: 300 K cooling: 10^ K cooling: 10^ K 
Top? 1.1x104 K Toe: 105 K Top? 1.5x104 K 


ps: 3.85x10?gcm? psp: ~7 x 108 gcm? (?) Pop: 3.x 107^ g cm? 


Although SF in the disk will also depend on the gas surface 
density profiles assumed, 
Kapferer sees much more SF in the tail! 


Harder questions: Magnetic fields 


e MEs in ICM, in galaxy . : Vollmer, Hmm -Masur 


| OBSERVATIONS di: k a».  LRP1 MODEL 


e Polarisation of galactic 
MFs on front edge 


DECLINATION (J2000) 


e Draping of ICM MEFs, 
protects tail from 
fragmentation 
(idealised initial 
conditions) 


Ruszkowski+ 
14 


Caution must be taken when posing 
questions and interpreting results! 


Generalizations of 
"How RPS affects galaxies" 
are dangerous!! 


Review of [theoretical] models 


= Annalisa Pillepich 
= Cosmological hydro sims 


The underlying numerical and physical ingredients 


Cosmological 
Gravity 

MHD 

Uniform Volume 
Sims for 


Galaxy Physics 


The model behind Illustris ING Annalisa Pillepich, Ringberg, Theory Review, 2017/12/11 


The underlying numerical and physical ingredients 


Cosmological 


The model behind Illustris ING 


Working Assumption: 
LCDM => Initial Conditions 


Components: 
Dark Matter ( {m articles) 
Dark Energy (global evolution) 


~250 Mpc (matter density projection) 


log,, (pdf) 


Pillepich et al. 2008 


Annalisa Pillepich, Ringberg, Theory Review, 2017/12/11 


The underlying numerical and physical ingredients 


Cosmological 


Gravity 


The model behind Illustris ING 


Working Assumption: 
LCDM => Initial Conditions 


Components: 
Dark Matter ( fm articles) 
Dark Energy (global evolution) 


Newtonian Equations 
in an expanding universe 
(non linear gravitational collapse) 


~250 comoving Mpc (matter density projection) 


Pillepich et al. 2008 


Annalisa Pillepich, Ringberg, Theory Review, 2017/12/11 


The underlying numerical and physical ingredients 


Hydrodynamics: 
+ Gas( Q} 


(H/He at the initial conditions) 


Euler Equations 


MHD 


The underlying numerical and physical ingredients 


Hydrodynamics: 
+ Gas ( ) 
(H/He at the initial conditions) 


Euler Equations 


MHD 


+ Seed Cosmic Magnetic Fields 
i.e. Maxwell Equations for perfect conductors 


z= 2.00 , 7 0.50 = 0.00 


10 Mpc 10 Mpc i 10 Mpc; 1 10 Mpa 


Pakmor et al. 2013, 2014, 2017 


Marinacti et al. 2016 Magnetic Field Strength, amplified from an initial 10^-14 Gauss to a few microGauss 


The model behind Illustris ING Annalisa Pillepich, Ringberg, Theory Review, 2017/12/11 


The underlying numerical and physical ingredients 


Feedback from stars 
+ STARS and BLACK HOLES => galactic outflows 


Formation and evolution of 
stars (SNla, SNII, AGBs) 

& their pollution of the 
inter-stellar medium 


(H, He, C, N, O, Ne, Mg, Si, and Fe) 


Pillepich et al. 2017 


Feedback from SMBHs 
=> galactic outflows & suppression of SF 


Galaxy Physics i ILEUS 


Sims for 


(via tables and including 
collisional excitation, collisional 
ionization, recombination, 
dielectric recombination and 
free—free emission) + UV 
background + metal line cooling 


The model behind Illustris ING 


Annalisa Pillepich, Ringberg, Theory Review, 2017/12/11 
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The underlying numerical and physical ingredients zc 


Uniform Volume 


The model behind Illustris ING 


e DARK MATTER 

resolution GAS - cells 

elements * STARS AREPO approach 
e BLACK HOLES 


Same resolution element mass across the whole box 
(e.g. DM particles mass = ~1046 Msun) 


Gas Mesh in a cosmological gas halo 


Full spatial and time resolution OLY OA T 
KEATING ane, 
adaptivity PL PA ae 


Thousands of galaxies and haloes! HIN 
j LR) 


O Q7 
Nelson et al. 2015 X 


Annalisa Pillepich, Ringberg, Theory Review, 2017/12/11 


The underlying numerical and physical ingredients 


Cosmological 
Gravity 

MHD 

Uniform Volume 
Sims for 


Galaxy Physics 


The model behind IllustrisTNG 


STELLAR LIGHT 


z=1.35 log, ,(M,)=11.5 SFR=162.6 sSFR-0.46Gyr 


THE HIERARCHICAL GROWTH OF GALAXIES, GALAXY MERGERS, 
COSMIC GAS ACCRETION INTO HALOES, TIDAL AND RAM 
PRESSURE STRIPPING, DYNAMICAL FRICTION etc ARE ALL 
"EMERGING" PROCESSES IN SIMULATIONS LIKE ILLUSTRIS/TNG 


GAS DENSITY 


Review of [theoretical] models 


= Dylan Nelson 
= ..inthe context of Illustris / TNG 


The original Illustris simulation: 


AREPO code (TreePM-Voronoi Hydro) 
100 Mpc cosmological volume 
Resolution: 1 kpc (109 M,,,, baryon) 
Galaxy mass range: > 10? Mau 

Halo mass range: < 10'^ M, 
‘Comprehensive physical model’ 


— Moster+ 2013 
—^ Behroozi4 2013 
== Kravtsov+ 2014 
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The original Illustris simulation: | « ses 
= AREPO code (TreePM+Voronoi Hydro) 

= 100 Mpc cosmological volume 

= Resolution: 1 kpc (109 M,,,, baryon) 
= Galaxy mass range: > 107? M,,, 

= Halo mass range: < 1014 Man 

= ‘Comprehensive physical model’ 
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disk galaxies 


Hir 


irregulars 


Takeaway differences: 
1. New BH feedback 
model: kinetic wind 
2. Revised galactic 
winds (Supernovae 
feedback) model 


Illustris TNG 


‘the next generation’ 


Volker Springel 
Lars Hernquist 
Annalisa Pillepich 
Ruediger Pakmor 
Dylan Nelson 


Rainer Weinberger 
Federico Marinacci 
Jill Naiman 


Shy Genel 
Paul Torrey 


ark Vogelsberger 


300 Mpc 


IllustrisTNG: are the results realistic? 
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IllustrisTNG: are the results realist 
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IllustrisTNG: are the results realistic? 


Galaxy color , 
bimodality 
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Galaxy evolution in Groups and Clusters: 
Prospects with TNG50 


Stellar Mass Function [ Mpc>3 dex! ] 
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Galaxy Stellar Mass [log Msun | 


Galaxy evolution in Groups and Clusters: 
Prospects with TNG50 


Stellar Mass Function [ Mpc-3 dex! ] 
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Galaxy Stellar Mass [log Msun ] 


Dylan Nelson, Ringberg, 11 Dec 2017 


Galaxy evolution in Groups and Clusters: 
Prospects with TNG50 


M32 SMC LMC 


Stellar Mass Function [ Mpc-3 dex! ] 
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Galaxy Stellar Mass [log Msun ] 


Dylan Nelson, Ringberg, 11 Dec 2017 
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Galaxy evolution in 


400 ckpe Groups and Clusters: 
Prospects with TNG50 
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Review of [theoretical] models 


= Yannick Bahé 
= ..inthe context of EAGLE 


The EAGLE simulations 


Many additional 50 and 25 cMpc simulations 

- Expore variations in supernova / BH 
feedback and other subgrid parameters. 

- Some groups in 50 cMpc volume 


Largest simulation: (100 cMpc)? 


- Resolution of ~2 x 10° Mo (baryons) 
- One "Reference" simulation model 
- Some (poor) clusters 


Supernova feedback in EAGLE 


log (tcooi / tcross, kernel) SN efficiency scalings 
(Dalla Vecchia & Schaye, 2012 (Crain et al., 2015) 
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Thermal stochastic feedback Feedback efficiency scaled with local quantities 
Energy is stored up until gas can be heated To overcome remaining numerical cooling, feedback 
to a temperature that is high enough for is made more efficient in high-density gas. Variation 


numerical cooling to be suppressed with metallicity models physical cooling losses. 


Supernova feedback in EAGLE 
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GSMF not sensitive to feedback scaling, but galaxy sizes are Crain et al., 2015 


EAGLE BH growth and feedback 


Modified Bondi accretion limit: 
Accretion disk: 


Maisc _ MBonditBondi 
viscous transport frre, Macer ~ = ~ EE UT ED 
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Free parameter Cvisc 


Describes unresolved structure of 
accretion disc (default Cvisc = 2rr). 


10019 May [Mc] 


Rosas-Guevara et al., 2015 


Higher Cvisc : lower viscosity, 
Schaye et al., 2015 


delayed BH growth 
Crain et al., 2015 
Black hole feedback is stochastic (like SN) 
BH feedback: ED 
But heating temperature is higher (default: AT = 
10% of accreted mass converted to energy, 1095 K) because gas density around black holes 
with assumed 15% coupling efficiency to gas 


is typically higher than in SF regions. 


Review of [theoretical] models 


= Current landscape -> collective discussion 


The Landscape of (Current) Cosmological (M)HD Simulations 


STNG (2017)... 
ILLUSTRIS (2014) 
EAGLE (2015). 


NIHAO (2015) 
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The current simulation landscape 
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The Landscape of (Current) Cosmological (M)HD Simulations 
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The current simulation landscape 
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The Landscape of (Current) Cosmological (M)HD Simulations 


TNG (2017)... 

ILLUSTRIS (2014) — 
EAGLE(2015)m. 5 
HORIZON-AGN (2014) | 
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MUFASA (2016) ^ 9 


MAGNETICUM (2017) ME 


Yonsei Clusters (2016) 
16 halos > 10114 Msun 


100 pkpc 


* Disk galaxies 


cc E 
30 haloes > 10414 Msun + 
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HYDRANGEA (2017) DMO Zooms 


:24 haloes » 10^14 Msun f . 
Adiabatic Runs 


Too-low res samples to study galaxies 


ZOOM IN SIMULATIONS of CLUSTERS 


It is hard to simulate very massive objects including the resolution and the physics 
ingredients that are needed to model also their member galaxies (and not just the central) 


The current simulation landscape Annalisa Pillepich, Ringberg, Theory Review, 2017/12/11 


The Landscape of (Current) Cosmological (M)HD Simulations 


1 + # Haloes (> M200c ) 


== TNG300 (m ~ 10’ Msun) 
7 TNG100 (m ~ 10° Msun) 
== C.EAGLE (m ~ 10° Msun) 
Rhapsody-G (m ~ 1084 Msun) 
Trieste (m ~ 108° Msun) 
= lh 


103 104 
Total Halo Mass [ M5006 ll Mn ] 


The current simulation landscape 


1013 


The mass regime of Fornax is not too 
much explored in zoom-in cosmological 
hydro simulations. 


In simulations like Illustris, EAGLE and 
the new TNG100, we have many tens of 
Fornax-like haloes and ~10 


In projects like C-EAGLE, there are 7 
haloes more massive than 10^15 Msun, 
so a few comparable to Coma. 


In TNG5O, we will have ~10 Fornax- 
like haloes with galaxies resolved 
with m res ~ a few 10^4 Msun 

and one Virgo mass-like object. 


Annalisa Pillepich, Ringberg, Theory Review, 2017/12/11 


Theory model questions (that we are asking ourselves, or that you should be asking us): 
e What are the strengths and weaknesses of current models in general? 
e For insights into galaxy evolution? For cluster galaxy members? 
* Are feedback models (i.e. AGN/ICM interaction) good enough, or too crude? 
* Are gas-dynamical models good enough, or too crude (i.e. ICM/ISM interaction)? 
e Can we do better: bridging the (res.) gap (i.e. zooms <-> boxes, idealized <-> SAMs) 
e Can models explain both ‘quiet gas stripping’ (gradual truncation) and ‘spectacular gas 
stripping’ (jellyfish, peri-center SF bursts) simultaneously? 
e What are the most important caveats in making obs. comparisons? Should we at all? 
e Are there un-considered obs. we should use to better constrain the models, or do the 
models have enough on their hands already? 
* Do group/cluster member constraints provide more than ‘field’ constraints? 
e Stepping back: do we understand DM sub-structure formation (& disruption!) processes 
enough to even model baryonic effects on top? 


Science questions (we) would like answers to: 
* How do different environmental effects impact satellite quenching (vs. mass/redshift)? 
e When (what distance) do satellites start to experience env. effects? 
e Does env. quenching depend on host halo mass, or not? (Just cosmic starvation?) 
* How does FB from central galaxies effect satellites? Directly, indirectly? 
e What direct evidence exists for the role of AGN FB? Distinguishable from SF FB? 
e What aspects of galaxies are affected by their environment? (e.g. AGN activity?) 
e Morphological evolution of satellites: important or not? 
* Metal enrichment of satellites: most important processes? Predictions? 


